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FOREWORD 
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Air  Force  Materials  Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio 
lt5l;33.  The  work  was  conducted  under  Project  73^2,  "Fundamental  Research 
on  Macromolecular  Materials  and  Lubrication  Phenomena",  Task  73^202, 
"Studies  on  the  Structure-Property  Relationships  of  Polymeric  Materials". 
The  program  was  administered  by  the  Air  Force  Materials  Laboratory. 
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ABSTRACT 

Failure  surface  of  laminated  composites  can  be  constructed  in  special 
strain  space  such  that  the  transformation  of  strain  components  becomes 
an  orthogonal  matrix.  This  construction  provides  a  convenient  means  of 
studying  strength  of  laminates  consisting  of  arbitrary  lamina  orientations. 
This  special  construction  of  failure  surface  can  be  based  on  the  maximum 
strain  theory,  the  tensor  polynomial  theory  or  other  failure  criteria  of 
the  lamina.  The  effect  of  nonlinearity  due  to  shear  on  the  failure 
surface  is  also  illustrated. 
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SECTION  I 
INTRODUCTION 

In  the  design  of  composites  it  is  frequently  necessary  to  know  the 
failure  properties  referred  to  a  reference  coordinate  system  other  than 
the  material  coordinates.  A  failure  surface  is  usually  constructed  in 
the  stress  or  strain  coordinates  which  are  the  components  referred  to 
the  material  symmetry  axes.  When  a  different  coordinate  system  is  used, 
the  resulting  failure  surface  cannot  be  obtained  from  the  original  one 
through  a  pure  rotation,  because  the  matrix  relating  the  stress  or  strain 
components  in  the  two  different  coordinates  is  not  orthogonal.  In  other 
words,  the  transformation  of  stress  or  strain  components  involves  a 
deformation  as  well  as  a  rotation,  which  follows  easily  from  the  Polar 
Decomposition  Theorem. 

Rotation  is  easy  to  visualize  and  so  amenable  to  graphical  representa¬ 
tion.  We  described  in  Reference  1  a  method  to  eliminate  the  deformation 
part  from  the  transformation  matrix.  What  makes  this  possible  is  the 
similarity  of  the  transformation  matrix  to  an  orthogonal  matrix.  Here  we 
apply  the  method  to  study  failure  surfaces  of  symmetric  laminated  com¬ 
posites  subjected  to  in-plane  loadings  only.  A  full  advantage  is  taken  of 
the  orthogonality  of  the  modified  transformation  matrix.  Further,  noting 
that  this  orthogonal  matrix  describes  a  rotation  about  the  axis  of  two 
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dimensional  dilatational  strain  p  in  the  modified  strain  coordinates  (or 
hydrostatic  stress  in  the  modified  stress  coordinates),  it  is  proposed  to 
represent  the  failure  surface  in  the  (ej  -  e2)/2-  e^ /2  plane  with  p  as  a 
parameter.  An  advantage  of  the  proposed  method  is  shown  via  its  appli¬ 
cation  to  analysis  of  initial  failure  of  a  laminate. 

Stress- strain  relation  of  most  unidirectional  laminae  can  be  assumed 
to  be  linear  except  in  shear.  A  comparison  is  given,  with  due  regard  to 
this  shear  nonlinearity,  between  the  stress  and  strain  criteria  of  failure 
when  they  both  are  described  by  polynomials  of  second  order. 
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SECTION  II 
METHOD 

It  was  suggested  in  Reference  1  to  use  \/~Z  o-fe  ^  in  place  of  the 
shear  stress  erg  because  this  makes  the  transformation  matrix  of  stress 
orthogonal.  The  same  result  was  shown  to  follow  for  strain  if  e^  is 
replaced  by  e^  /V  2  .  A  new  set  of  coordinates  (o^i*,  ^2*,  was 

introduced  in  which  the  transformation  describes  a  rotation  about 
axis  through  an  angle  -20,  where  0  is  the  angle  of  rotation  of  coordinates 
(Xi»  X2)  Figure  1, 

The  modified  strain  components  e^^  are  related  to  the  original  com¬ 
ponents  e^  by 


ei 

(61  +  62)/"'/  2 , 

62*  = 

-(ei  -  ez)l\  2, 

66*  = 

To  make  terms  on  the  right  side  more  familiar  we  divide  them  by 
to  obtain 

e°  =  (p,  q,  r)T  =  He  (2) 

Contracted  notation  is  used. 
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where 


Note  that  the  coordinate  system  (p,  q,  r)  is  left-handed. 

Quantities  p,  q,  r  given  by 

1 

p  =  (ei  +  e2)/2,  q  =  (ej-  &z)l2,  r  =  e^/Z,  (U) 

are  all  familiar  if  we  consider  the  Mohr's  circle  representation  of  a 
strain  state,  Figure  2.  Moreover,  we  can  specify  (e^,  e^)  and  (p,  q)  in  the 
same  coordinate  plane  if  we  change  the  scale  of  ei  and  eg  by 

©1  =  ©1  /V^  2,  62  =  62  /\^  2^ 

The  relation  is  shown  graphically  in  Figure  3.  Transformation  of  the 
modified  coordinates  (p,  q,  r)  follows  from  that  of  (ei,  62,  -©e): 


where  primed  quantities  are  referred  to  new  coordinates  (yj ,  X2)  • 
Equation  (6)  is  nothing  but  a  simple  coordinate  transformation  showing 
that  (p',  q',  r')  is  related  to  (p,  q,  r)  through  a  rotation  of  angle  20  about 
p  axis.  Figure  it. 
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Figiire  2.  p,  q,,  r  in  Mohr's  Circle 
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Figure  3-  Relation  Between  (e^,  e^)  and  (p,  q) 
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Failure  is  characterized  by  a  closed  surface  in  the  stress  or  strain 
coordinates.  Transformation  between  e^'  and  e^  for  example,  includes 
a  deformation  as  well  as  a  rotation.  Therefore,  the  shape  of  failure 
surface  is  altered  by  this  transformation,  Figure  5.  However,  when  the 
failure  surface  is  drawn  in  the  (p,  q,  r)  coordinates,  its  shape  does  not 
change  because  the  coordinate  transformation  results  only  in  a  rotation 
about  the  p  axis. 

Mathematically,  a  failure  surface  may  be  described  by  a  function  g 
such  that 

g(e)  =  1  (y) 

The  equation  for  the  same  surface  in  the  (p,  q,  r)  coordinates  is  then 
given  by  a  new  function  g°  which  is  related  to  g  by 

g°(eO)  =  g(H-ieO)  =  1  (8) 

Equation  (8)  is  useful  in  obtaining  the  new  function  g°,  especially  when 
the  surface  is  described  by  a  polynomial.  It  should  be  noted,  however, 
that  the  method  is  not  restricted  to  such  cases  only.  A  simple  graphical 
method  is  possible  when  the  maximum  strain  criterion  is  used.  Figure  6 
shows  such  graphical  construction.  It  also  illustrates  the  way  the  projec¬ 
tion  on  the  q-r  plane  is  obtained.  This  will  be  explained  in  detail  in  the 
following  section. 

The  same  reduction  follows  for  a  failure  surface  in  the  stress  space  if 
the  stress  vector  a  =  (eri  ,  era,  crfeF  is  replaced  by 
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0-°  =  (p,  q,  r)*^  =  Ko- 

where 
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(9) 


(10) 


Our  discussion  is  hereafter  restricted  to  strain  criteria  because  our  main 
objective  is  to  apply  the  method  to  laminated  composites  in  the  framework 
of  classical  plate  theory. 
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SECTION  III 
APPLICATIONS 

1.  UNIDIRECTIONAL  LAMINAE 

There  are  various  criteria  for  failure  of  composites.  These  theories 

are  summarized  in  Reference  2.  Although  the  present  method  is  not  restricted 
to  any  particular  criteria,  ve  restrict  our  discussion  to  the  maximum 

strain  criterion  and  the  second  order  tensor  polynomial  criterion. 

In  the  maximum  strain  criterion,  the  allowable  range  of  strains  are 

given  by 

El'*’  >  ei  >  El",  Ea"*"  >  ea  >  Ea", 

,  _  (11) 

E6'^  >  eg  >  Eg 

In  the  modified  coordinates  (ij,  ea,  r)  the  corresponding  limit  values 
are 

Ei+/V^  >  e"!  >  Ei'/vX  Ea+ZV  2  >  ez  >  Ea’/'^Z, 

Ee'^/Z  >  >  Efi'/Z  (12) 

Equation  (ll)  implies  the  absence  of  interaction  among  limiting  strain 
components.  Therefore,  it  is  possible  to  obtain  those  limiting  values 
from  the  limiting  stress  components  determined  by  load- controlled  tests 
if  pertinent  constitutive  equation  is  known.  Most  composites  are  almost 
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linear  up  to  failure  in  longitudinal  and  transverse  tension  and  compres¬ 
sion.  Thus,  it  follows  that 


Ei+  =  Xi+/El, 
E2+  =  X2+/E^, 


Ef  =  xr/E^, 
E2‘  =  X2'/E^, 


(13) 


where  Xi  and  X2  are  the  failure  stresses  and  Ej^  and  Erp  are  the  longi¬ 
tudinal  and  transverse  Young's  moduli,  respectively  (Reference  3). 

The  difficulty  in  relating  the  limiting  strain  components  to  the  limiting 
stress  components  lies  in  shear,  because  of  highly  nonlinear  behavior. 

The  procedure  we  follow  here  is  derived  from  current  practice  which  chooses 
a  fictitious  maximum  strain  assuming  a  linear  behavior.  This  may  be 

plausible  if  we  note  that,  the  classical  theory  of  laminated  composites  is  based 
on  the  linear  stress-strain  relation.  A  recent  study  (References  k,  5)  showed 

that  the  effect  of  shear  nonlinearity  on  the  behavior  of  laminates  is  negligible 
under  uniaxial  tension  as  far  as  laminates  contain  layers  of  0°  orientation. 
The  nonlinearity  appears  most  in  ±45°  laminate  as  experimentally 
observed  and  theoretically  proved.  In  any  event  we  here  choose  the 
truncated  value  as  the  limiting  shear  strain  so  that 

E6+  =  -E6-  =  Xfi/G  (14) 

where,  again,  Xg  is  the  failure  stress  and  G  is  the  shear  modulus.  The 
first  equality  which  states  the  identity  of  positive  and  negative  shears 
follows  from  the  material  symmetry  of  unidirectional  laminae. 
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Failure  stresses  of  unidirectional  graphite/epoxy  and  boron/epoxy 
composites  are  given  in  Reference  6  .  The  corresponding  strains  are 
determined  using  Equations  13  and  l4,  and  are  listed  in  Table  I.  Graphical 
representation  of  failure  surfaces  is  given  in  Figures  6  through  9- 

The  procedure  to  be  followed  is  described  below: 

a.  Obtain  limiting  values  of  ei ,  ez,  eg  directly  from  experiment 
or  from  limiting  stress  components  through  Equations  13  and  l4. 

b.  Calculate  limiting  values  of  ej ,  ez,  r  from  Equation  12. 

c.  Draw  the  failure  surface  in  the  'ej -"ea  plane. 

d.  Draw  p  and  q  axes,  p  axis  bisects  the  angle  <  ejoea  and  q 
axis  the  angle  <"610  (-"ea). 

e.  To  find  the  failure  surface  in  the  q-r  plane  corresponding  to  a 
particular  value  p  =  c,  draw  a  straight  line  p  =  c.  The  points  at  which 
this  line  intersects  the  failure  surface  in  the  p-q  plane  give  the  limiting 
values  of  q. 

f.  The  corresponding  failure  surface  is  then  given  as  a  rectangle  whose 
vertical  sides  pass  through  those  determined  points  on  the  q  axis  and 
whose  horizontal  sides  are  given  by  r  =  Eg'^/2  and  r  =  Eg  /2. 

Another  way  of  characterizing  the  failure  of  unidirectional  composites 
is  provided  by  a  second  order  tensor  polynomial  of  the  form  (Reference  7)* 


*  Summation  is  implied  over  any  repeated  indices. 
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TABLE  I 

LIMIT  VALUES  OF  STRAINS 


Failure 

strain 

Material 

X  10^ 

-El' 

X  10^ 

E2+ 

X  10^ 

-Ea" 

X  10^ 

E6 

X  10^ 

High-Strength 

Graphite  /  Epoxy 

V£  =  0.6 

8.57 

8.57 

4.71 

17.65 

18.46 

Intermediate-  Strength 
Graphite  /  Epoxy 

V£  =  0.6 

9.41 

9.41 

4.41 

14.71 

15.38 

High- Modulus 

Graphite  /  Epoxy 

V£  =  0.6 

4.40 

4.00 

2.35 

11.76 

13.85 

Boron/ Epoxy 

V£  =  0.5 

6.40 

11.77 

3.85 

14.81 

21.86 
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Figure  6.  Maximum  Strain  Failure  Siarface  for  High-Strength 
Graphite/Epoxy 
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Figure  7-  Maximum  Strain  Failiore  S\rrface  for  Intermediate- 
Strength  Graphite /Epoxy 
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Figure  8.  Maximum  Strain  Failure  Surface  for  High-Modulus 
Graphite/Epoxy 
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g(e)  =  Gie-  +  Gjjeiej  -  1 


In  the  modified  strain  coordinates  the  same  failure  surface  is  described 


by  a  new  function  g° 


gO(eO)  .  Gi°e,o  +  G.j°e.°e.o  .  ^ 


where 


Gi°  =  h:/  g  , 


G-.°  =  H‘!  G,  „ 
ij  ki  ij  kf 


Using  Equation  (3)  for  the  matrix  in  Equation  (17)  results  in 

^  J 


Gi  +  G; 


Gl  -  Ga  >  , 


Gii  +  2Gi2  +  G22 


symmetric 


Gii  -  G22 


2(Gi6  +  G26) 


Gji  -  2Gi2  +  G22  2(Gi6  ■  G26) 


For  unidirectional  laminae  the  material  symmetry  requires  that 


Ge  -  Gi6  -  G26  -  0 


Therefore,  Equation  (16)  reduces  to 
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Since  a  failure  surface  should  be  closed  from  physical  considerations, 
Equation  (20)  should  describe  an  ellipse  and  so 

a^>0,  b^>0,  ®>kp>0  (22) 

The  ratio  of  major  axis  to  minor  axis  remains  constant  whereas  the  center 
moves  along  q  axis  depending  on  the  value  of  p.  The  lengths  of  axes  also 
depend  on  p.  In  terms  of  Equation  (22)  yields  the  well  known  conditions 

Gii°G230  >  Gl2°^  G22°G66°  >  0  (23) 

and 

Gi°%Gii°G2°Vg22°  +  4G„°-  2Gi2°(Gi°G2°+  2Gi2°)/G22°  >  0  (2^+) 

The  bo\mds  for  p  then  follow  as 

<  P  <  Pma*,  <25) 
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p  .  =  (B  -  C)/(2A), 

^min  '  !  \  ’ 


p  =  -  (B  +  C)/(2A) 
^max  '  ’  ^  ' 


(25) 

Cont'd. 


where 


A  =  Gi2°  /G22O  -  G„° 


B  =  [Gi°%  GiiOG2°Vg22°  +  4Gii°-  2Gi2°(Gi°G2°  +  2Gi2°)/G22°] 


G  =  G2°Gi2°/G22°  -  Gi° 

The  strength  tensors  G^,  Gj^j  in  the  strain  criterion  can  be  related  to 
the  strength  tensors  Fj^j  in  the  stress  criterion  if  the  material  is 

linear  elastic  (Appendix  I): 


Gi  =  Q.jFj, 


"^ij  ^ki^fj^ki 


where  Q^j  is  the  reduced  stiffness  tensor.  The  above  equation  is  necessary 
because  most  available  data  are  for  Fj^  and  Fjj. 

As  an  example  we  take  the  failm’e  data  reported  in  Reference  8  for 
glass/epoxy  scotch-ply  1002  lamina: 


Fi  = 

-0.003  (ksi)-^ 

F2  = 

0.295  (ksi)-' , 

Fn  = 

0.063  X  10“^  (ksi)-^. 

II 

17. 15  X  10"^  (ksi)- 

F66  = 

10.85  X  10-^  (ksi)-^. 

Fi2  = 

0.  15  X  10-^  (ksi)- 

The  stiffness  matrix  of  the  same  material  is 


5.004  0.083  0 

[Q^j]  =  0.083  1.668  0 

0  0  0.704 


X  10^  ksi 
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From  Equation  (27)  follow  G.  and  G.  : 

1  ij 


Gi  =  9.475 
G„  =  1.82  X  10^ 

G66  =  5. 38  X  10^ 

Finally,  substituting  Equation  (30) 


II 

O 

510.29 

II 

r— 1 

o 

56.88  X 

10 

G66°  = 

21.51  X 

10 

Ga  = 

49 1 . 8 1 

G22  = 

47.76  X  10^ 

(30) 

Gi2  = 

3. 65  X  10^ 

into  Equation  (18)  yields 

G2°  = 

-482.34 

G22°  = 

42.27  X  10^ 

(31) 

Gi2°  = 

-45.94  X  10^ 

The  variables  characterizing  the  failure  surface  can  now  be  obtained 
by  substituting  the  above  data  into  Equation  (21): 

a  =  4.864  X  lO'^kp,  b  =  6.818  x  IQ-^kp, 

=  1.087  p  +  5.705  x  lO'S 

^  (32) 

kp  =  (-6.951  X  10^  +  22.928p  +  2.376)^/^. 

Pmin  =  ^  10-^  p^^^  =  20.21  x  lO'^ 


The  failure  surface  is  drawn  on  the  q-r  plane.  Figure  10.  Each  ellipse 
represents  Equation  (20)  for  a  fixed  dilatational  strain  p. 

It  goes  without  saying  that,  whereas  the  graphical  method  is  a  more 
viable  means  of  constructing  failure  surface  in  the  case  of  maximum 
strain  (or  stress)  criterion,  one  has  to  resort  to  the  equation  of  failure 
in  the  case  of  polynomial  criterion  (Equation  15). 
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Figure  10.  Polynomial  Failure  Surface  for  Glass/Epoxy  Scotch-Ply  1002 
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Now  that  one  has  the  appropriate  failure  surface  drawn  in  the  q-r 
plane,  any  off-axis  failure  can  easily  be  checked.  This  will  be  illustrated 
in  the  following  paragraph. 

2.  LAMINATED  COMPOSITES 

According  to  the  classical  plate  theory  for  thin  laminates,  the  in- plane 
strains  are  assumed  to  be  the  same  in  all  constituent  layers.  This  makes 
it  possible  to  draw  failure  surface  for  each  layer  in  the  same  strain 
coordinates  without  modification.  Once  the  state  of  critical  strain  is 
known,  the  state  of  stress  in  each  layer  can  be  calculated  by  using  appro¬ 
priate  constitutive  equations,  e.g.,  (r^  =  Q^jCj. 

The  above  argument  applies  at  most  to  the  initial  failure  of  laminate, 
i.e.,  the  failure  of  the  weakest  layer  in  the  laminate.  The  behavior  after 
the  initial  failure  is  not  well  understood  at  this  time.  One  of  the  approaches 

to  this  problem  is  to  remove  the  layers  that  have  failed  from  the  laminate 
by  setting  the  stiffness  of  the  layers  equal  to  zero  (Reference  9)j  the 

other  assumes  stress-strain  relations  with  negative  slope  for  the  failed 

layers  (Reference  10 ).  In  any  event,  understanding  of  the  initial  failure 

is  important  because  it  forms  one  of  the  bases  of  current  design  criteria 

(References  6,  11 ). 

One  more  point  to  be  noted  is  that  we  do  not  consider  any  effects  of 

interlaminar  stresses.  The  interlaminar  stresses  become  significant  along 
free  edges  because  there  they  have  large  magnitude,  depending  on  the  stacking 
sequence  or  lay-up  of  the  laminae  (References  12,  13).  However,  a  few 
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thicknesses  of  the  laminate  away  from  the  free  edges  the  classical  plate 
theory  provides  a  reasonable  solution. 

To  he  precise  we  formulate  the  criterion  for  the  initial  failure  as 
follows  (Reference  li+):  Let  g  be  the  failure  function  of  a-th  layer. 

Then  the  initial  failure  of  laminate  occurs  when  the  strain  e  satisfies 
the  condition 

n:iaxJga(e)l  =1  (33) 

a  ^  ~  ^ 

Graphically,  this  is  the  equation  for  the  inner  envelope  of  failure  surfaces 
of  constituent  layers. 

As  we  have  seen,  the  failure  surface  of,  e.g.,  (j)°  layer  is  obtained 
from  the  failure  surface  of  0°  layer  by  keeping  the  surface  fixed  while 
rotating  q,  r  axes  through  angle  -2<j).  Note  that  the  layer  orientation  angle 
cj)  is  equal  to  -0,  because  9  is  measured  away  from  the  layer  symmetry 
axes  (Figure  6). 

As  an  illustration  consider  a  high-strength  graphite /epoxy  laminate  of 
[0°/30°]g  lay-up  configuration.  The  lamina  failure  surface  is  shown 

for  the  0°  layer  in  Figure  6.  Suppose  the  state  of  strain  resulting  from 
stress  analysis  is 

e^^  =  4x10“^,  By.  =  -4x10"^,  e^y  =  12x10"^  ^^4) 

The  corresponding  p,  q,  r  follow  from  Equation  4  by  using  the  above 
values : 
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p  =  (e^  +  e^)/2  =0,  q  =  (e^  -  e^)l2  =  4  x  10"^, 

r  =  exy/2  =  6  X  10-^  (35) 

The  state  of  strain  is  denoted  by  the  point  Q  for  the  0°  layer.  Since  we 
use  the  same  failure  surface,  the  same  state  of  strain  is  given,  for  the  30° 
layer,  by  the  point  Q'  that  is  obtained  by  rotating  0  Q  through  angle  -60°. 
The  failure  surface  corresponding  to  p  =  0  is  shown  by  dashed  lines. 

Since  both  Q  and  Q'  are  inside  the  region  enclosed  by  those  dashed  lines, 
both  layers  will  not  fail.  It  is  easily  seen  that  under  the  same  state  of 
strain,  the  range  of  cjj  without  failure  is 

-74°  <  24)  <  186°,  i.e.,  -37°  <  4)  <  93°  (36) 

One  of  the  advantages  of  using  (p,  q,  r)  coordinates  in  conjunction  with 
failure  is  that  one  constructs  the  failure  surface  only  once  no  matter  how 
many  layers  a  laminate  consists  of  insofar  as  all  the  layers  are  made  of 
the  same  material. 

A  disadvantage  of  working  with  failure  surfaces  in  the  stress  coordinates 
stems  from  the  fact  that  the  lamina  stresses  vary  from  layer  to  layer  even 
when  there  are  no  bending  effects.  Thus  the  relation  between  the  lamina 
stresses  and  the  laminate  stresses  must  be  known  in  order  to  express  any 
failure  criterion  in  terms  of  the  laminate  stresses. 

In  the  laminate  coordinates,  the  laminate  stress- strain  relation  is 
given  by  Equation  3T  (Reference  15)- 
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[An]  e. 


where  [-A-ij]  is  the  laminate  stiffness  matrix.  Noting  that  the  strain  com¬ 
ponents  in  the  lamina  coordinates  are  related  to  those  in  the  laminate 
coordinates  by 


j 

cos  ^cj) 

sin^cj) 

(sin 2  <|))/2 

[^fjl  'y  • 

II 

sin^<j) 

cos 

-(sin2  4>)/2 

■*xy) 

-  s  in  ^  cj) 

sin^cj> 

cos  2  cj) 

we  can  obtain  the  relation  between  the  lamina  stresses  and  the  laminate 
stresses: 


tA,]  N, 


=  °ik  A-; 


where  i®  inverse  matrix  of 

To  compare  the  failure  behavior  of  a  lamina  as  an  independent  layer 
and  as  a  constituent  layer  in  a  laminate,  it  is  necessary  to  use  the  same 
reference  coordinates.  Thus  we  need  the  equation 


[Tij]  by  .  [Trj) 


2  , 

COS  9 

sin^cj) 

sin  2  cj) 

sin^ej) 

2  , 

cos  9 

-  sin  2  cj) 

(sin2  4>)  /2 

(sinZ  c())/2 

cos  2  cj) 
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for  the  -30°  layer. 

For  the  purpose  of  comparison,  we  choose  the  maximum  stress 
criterion  formulated  by 


> 

O'! 

>  xT  , 

X2  ^0'2>X2  , 

X6 

> 

0-6 

>  -X6 

{h6) 

For  the  given  boron/epoxy  composite,  the  limit  values  are 


II 

192  ksi. 

xf 

=  -353  ksi. 

X2^  = 

10.4  ksi. 

X2' 

=  -40  ksi. 

X6  = 

15.3  ksi 

As  an  independent  layer,  failure  surfaces  are  shown  in  Figure  11  for  0°, 

+30°,  and  -30°  orientations.  Since  =  0  in  the  case  under  consideration, 

xy 

and  since  =  -X^"  =  X^  ,  the  failure  surface  for  +30°  orientation  is 

the  same  as  that  for  -30°  orientation. 

Moreover,  for  ±30°  orientations  the  possible  failure  mode  is  either 
transverse  or  shear  failure.  In  the  figures  the  first  letter  stands  for 
longitudinal  or  transverse  and  the  second  for  tension  or  compression. 

S  denotes  shear. 

Failure  behavior  of  the  same  lamina  as  a  constituent  layer  in  the 
[0°/±  30° ]g  laminate  is  quite  different,  as  shown  in  Figure  12.  Again,  +30° 
and  -30°  layers  behave  identically,  because  =  0.  Note  that  even  the 

0°  layer  has  a  failure  surface  which  is  rotated  as  well  as  deformed  from 
the  original  shape. 
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Failure  Surfaces  of  Boron/Epoxy  Lamina  as  Independent  Layers 
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Figiire  12.  Failure  Surfaces  of  Boron/Epoxy  Lamina  as  Constituent  Layers  in  [0°/±30°]  Laminate 
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SECTION  IV 

COMPARISON  BETWEEN  STRESS  AND 
STRAIN  FAILURE  CRITERIA 

In  the  preceding  section  limiting  strains  were  obtained  from  failure 
stresses  under  the  assumption  of  linear  stress- strain  relation.  It  was 
noted  also  that  those  are  not  the  actual  values  but  fictitious  ones  because, 
frequently,  the  relation  is  not  quite  linear.  Unidirectional  composites  are 
observed  to  behave  nonlinearly  ,  especially  in  longitudinal  shear.  Recently, 
an  attempt  has  been  made  to  describe  analytically  this  shear  nonlinearity. 

A  reasonably  good  agreement  was  shown  between  the  analytical  predictions 
and  experimental  data  when  a  term  of  fourth  order  in  shear  stress  was 
added  to  the  plane-stress  complementary  energy  of  the  linear  elasticity. 

Physically,  failure  criteria  in  stress  and  in  strain  should  be  convertible 
to  each  other  because  they  express  the  same  phenomena.  What  prevents 
this  is  the  lack  of  both  exact  criteria  and  accurate  stress- strain  relation. 

It  is  shown  in  Reference  7  that  a  stress  tensor  polynomial  of  second  order 
provides  a  reasonably  accurate  failure  criterion.  Then  the  question  that 
might  be  asked  is,  how  does  it  compare  with  a  strain  tensor  polynomial 
of  the  same  order?  An  answer  to  this  question  will  be  attempted  by  taking 
a  stress-strain  relation  of  the  form  suggested  in  Reference  1+  . 

The  stress- strain  relation  we  shall  use  is  written  as 
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^1  1 

{  O'! 

^2  ) 

-  [Sij]  0-2 

66  1 

1 

(48) 


0'6 


where  the  linear  compliance  matrix  [%j]  is  given  by 

S, 


["ij]  = 


Sji  Si2  0 

Si  2  S22  0 

0  0  Sgfe 


(49) 


Then,  substituting  Equation  (48)  into  Equation  (l5)  yields 
f  (a-j_) 


Fieri  +  F2cr2  +  FucTj  +  ZFizcricrz  +  F 22  CTg 


+  F(,(,  (I  +  a.  (T^)^  (T^  =  1 


(50) 


where  F^  and  F^j  are  defined  by 


F.  =  S..G.,  F..  =  i=1.2 


and 


n  S5666 


(51) 


(52) 


Now  that  the  strain  polynomial  is  converted  to  a  stress  polynomial, 
the  stage  is  set  to  compare  the  stress  and  strain  criteria  of  failure  when 
they  both  are  described  by  polynomials  of  second  order.  For  this  purpose 
the  stress  polynomial  given  in  Appendix  I  is  rewritten  as 
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=  Fjo-i  +  F20-2  +  FnCi  +  2Fi2  ctictz  +  F 22  cr2  +  Fg^  tr^  =  1  (53) 

Five  coefficients  except  F^^  are  determined  from  five  different  tests 
(guiding  experiments  (Reference  16)  in  the  Oi  —  O2  plane:  tension  and 
compression  in  the  longitudinal  and  transverse  directions  and  a  biaxial 
loading).  Therefore,  it  follows  that 


F.  =  F.,  F..  =  F..; 

1  1  ij  ij’ 


i,  j  =  1,2 


(5it) 


However,  in  order  for  Equations  (50)  and  (53)  to  have  the  same  intercept 
±  Xg  along  the  erg  axis,  the  following  relation  should  hold: 


^66=  ^66 
(1  +  aXfe")" 


(55) 


Thus  the  failure  surface  in  the  stress  coordinates  representing  the  strain 
criterion  of  the  form  (Equation  15)  can  be  written  as 


f((r^)  =  Fj  o-j  +  F20-2  +  FiiOj^  +  2Fi2cri(r2  +  F zz 


(56) 


+  F 


66 


1  +  a  erg 
,  1  +  aXg  ' 


-6  =  1 


Suppose  that  (tri ,  <tz>  (rg)  and  (eri,  0-2,  erg)  satisfy  Equations  (53) 
and  (56),  respectively.  Then  we  have 


O-g  I  1  +  Q  O-g^ 


0'6 


1  +  aXg 


(57) 


The  above  equation  implies  that 
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_ 

> 

— 

< 

< 

0'6 

as 

<^6 

(58) 


In  other  words,  the  strain  criterion  predicts  a  higher  failure  stress  of 
shear  in  combined  loadings  when  the  applied  shear  stress  is  less  than  the 
failure  stress  in  pure  shear,  and  vice  versa. 

Figures  13  and  ik  show  ellipses  described  by  Equations  (53)  and  (56). 
The  data  except  a  are  for  graphite  fiber  (Morganite  II)  reinforced  epoxy 
composite  lamina  fabricated  by  Whittaker  Corporation: 


Fi  =  -0.003  (ksi)-^  Fa  =  0.105  (ksi)-^ 


6  - 


0, 


Fii  =  0.065  X  10-Mksi)-^,  F22  =  8.72  X  10-Mksi)-% 


F66  =  9.07  X  10-^  (ksi)-^ 


Fi2  =  0.2  X  10-^  (ksi)- 


The  values  of  a  are  taken  from  Reference  t  only  for  illustration  purposes. 
However,  the  actual  value  for  the  chosen  composite  is  expected  to  fall  within 
the  range  considered.  It  seems  that  with  higher  values  of  a  the  strain 
failure  criterion  tends  toward  the  maximum  stress  criterion.  At  present 
there  is  no  conclusive  experimental  evidence  which  can  support  the  choice 
of  one  over  the  other.  In  any  case,  the  difference  seems  rather  small. 
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SECTION  V 
SUM4AEY 

Current  practice  of  constructing  the  initial  failure  surface  of  laminates 
is  rather  complicated  and  usually  needs  the  help  of  a  computer.  More¬ 
over,  in  the  design  of  laminated  composites  it  is  frequently  necessary  to 
change  the  lay-up  configuration.  Thus,  for  each  case  a  failure  surface 
must  be  determined  anew. 

In  the  proposed  method  coupled  with  the  Kirchhoff  assumption,  the 
failure  surface  is  constructed  only  once  for  a  layer  with  given  direction, 
usually  0°  orientation.  Failure  surfaces  for  the  other  off-axis  layers  are 
obtained  from  that  for  the  0°  layer  simply  through  an  appropriate  rotation 
about  the  axis  of  dilatational  strain.  Alternately,  we  may  hold  the  failure 
surface, fixed  while  rotating  the  strain  coordinates  through  a  pertinent  angle. 

The  method  is  based  on  the  observation  that  the  transformation  matrix 
for  the  strain  components  becomes  orthogonal  when  the  shear  strain  eg 
is  replaced  by  eg /V^TT  For  the  stress  one  has  to  multiply  the  shear  stress 
by  \''z.  Consequently,  it  is  proposed  to  use 

p=-|-(ei  +  e2),  q=^(ei-e2),  ^  ~  ~z' 

in  place  of  ei ,  Bz,  eg  .  The  above  are  the  quantities  we  are  familiar  with 
in  the  Mohr's  circle  representation  of  a  strain  state.  Then  it  is  easily 
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verified  that  there  exists  a  one-to-one  correspondence  between  the 
rotation  of  (xi »  X2)  ^.nd  that  of  (q,  r) . 

One  of  the  disadvantages  may  be  that  uniaxial  strengths  are  not 
apparent  in  the  proposed  method.  This  can,  however,  be  resolved  by 
studying  the  relation  between  p  and  q,  and  if  failure  surfaces  are  available 
in  the  q-r  plane  for  a  reasonable  number  of  p  values. 

Using  the  constitutive  equation  suggested  in  Equation  1+  to  describe  the 
nonlinearity  in  shear,  a  comparison  was  made  between  the  stress  and  strain 
criteria  of  failure  when  they  both  are  described  by  polynomials  of  second 
order.  Although  there  is  a  noticeable  difference  in  combined  stress  state, 
no  experimental  data  are  available  yet  to  favor  one  over  the  other.  Con¬ 
sidering  both  the  difficulty  of  working  with  nonlinear  equations  and 
experimental  scatter,  the  difference  may  be  neglected. 
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APPENDIX  I 

MODIFIED  STRENGTH  TENSOR  COMPONENTS 


Strength  failure  criterion: 


f  (o-) 


FiUi 


F- .  cr-  cr. 
ij  1  J 


1. 


Modified  criterion: 


where 


f°(o-°)  =  F°o-°  +  F.°o-°a-°  =  1 

11  ij  1  1 


Fji  +  2F  12  +  F22 


Fn  ■  F  22 


Fii  -  ZFi2  +  F 


22 


symmetric 
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APPENDIX  II 

USE  OF  FAILURE  SURFACES  IN  q-r  PLANE 

Suppose  the  elastic  stress  analysis  of  a  laminated  composite  has 
yielded  the  following  strain  components  referred  to  the  laminate  reference 
axes : 

( e  ,  e  ,  e  ) . 

'  x’  Y  xy' 

Failure  surfaces  of  0°  layer  corresponding  to  various  values  of  p  are 
available  in  the  q-r  plane.  Failure  of  each  layer  under  the  given  state  of 
strain  can  then  be  checked  by  following  the  procedure  described  below 
(Figure  6) : 

a.  Calculate  p,  q,  r  using  Equation  (i;),  i.e., 
p  =  (e^  +ey)/2,  q  =  (e^-  ey)/2,  r= 

b.  Choose  the  failure  surface  corresponding  to  the  calculated  value 
of  p  in  the  q-r  plane. 

c.  Find  the  point  Q  with  coordinates  (q,  r)  and  draw  a  straight  line 
OQ  connecting  the  origin  O  and  the  point  Q. 

d.  For  a  cj)°  layer,  find  the  point  Q'  such  that  OQ'  is  obtained  by 
rotating  OQ  clockwise  through  angle  2(j). 

e.  If  Q'  is  inside  the  region  enclosed  by  the  chosen  failure  surface, 
the  layer  is  safe.  Otherwise,  the  layer  has  failed. 
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